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Abstract. We investigate the mass functions of different elements of the Large Scale Structure - walls, pancakes, 
filaments and clouds - and the impact of transverse motions - expansion and/or compression - on their statistical 
characteristics. Using the Zel'dovich theory of gravitational instability we show that the mass functions of all 
structure elements are approximately the same and the mass of all elements is found to be concentrated near 
the corresponding mean mass. At high redshifts, both the mass function and the mean mass of formed elements 
depend upon the small scale part of the initial power spectrum and, in particular, upon the mass of dominant 
fraction of dark matter (DM) particles. These results generalize the Press-Schechter approach and are used to 
obtain independent estimates of probable redshifts of the reionization and reheating periods of the Universe. We 
show that the transverse motions do not significantly change the redshift evolution of the observed mass function 
and the mean linear number density of low mass pancakes related to absorption lines in the spectra of the farthest 
quasars. We apply this approach to the observed Lyman-a clouds and obtain direct estimates of the variance 
of initial density perturbations and the shape of initial power spectrum on small scale. In turn, these estimates 
restrict the mass of dominant fraction of DM particles. 

Key words, cosmology: large-scale structure of the Universe — dark matter - galaxies: formation. 
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J> \ 1. Introduction However, this relation is not yet clear and a quantitative 

■ description of the LSS elements only recently got under 

■ In recent years numerical simulations are being more of- way Now the analysis of both simu l a ted and observed 
g . ten used to investigate the process of formation and evo- catalogues of ga l axi es is focused on the discussion of the 

lution of the Large Scale Structure (LSS) of the Universe. corre l a tion function, £(r), and the power spectrum, p(Jfe), 

This trend is strongly stimulated by the rapid growth of wMe other characteristics of the LSS are not determined, 

computer facilities and it provides better and better high Thig limited description demonstrates that numerical sim- 

quality simulations (see e.g., Benson et al. 2001; Smith ulations cannot subs titute theoretical models of structure 

et al. 2002; Frenk 2002). Modern simulations performed formation which provide the basis for much more detailed 

in large boxes and with high resolution can take into ac- quantitative description of both observed and simulated 

count simultaneous action of many important physical fac- j ar g e scale matter distribution 
tors and they substantially increase available information 

about the process of LSS formation and evolution. Theoretical models reveal the main factors that in- 
At small redshifts the LSS is observed in the large fluence the process of structure formation and evolution 
galaxy surveys as a system of filaments and walls. At and clarify links between them and the measured char- 
high redshifts the LSS is observed mainly as the Lyman- acteristics of the LSS. This opened up a possibility of 
a clouds identified through absorption lines in spectra of formulating new and promising approaches of statistical 
the farthest quasars. All these elements of the LSS are description of both simulated and observed characteris- 
quite well reproduced in numerical simulations what indi- tics of the LSS such as, for example, a set of mass func- 
cates their close relation with the initial power spectrum, tions of structure elements, their separations and so on. 
These characteristics can be measured with such powerful 
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(Barrow, Bhavsar feSonoda 1985, van de Weygaert 1991; 
Demiariski et al. 2000 hereafter DDMT; Doroshkevich et 
al. 2001) and the Minkowski Functional method (see, e.g., 
Kerscher 2000) among others. The well known example 
of this approach is provided by the Press-Schechter rela- 
tion for the mass function of high density clouds which has 
been widely used and discussed during the last thirty years 
(see, e.g., Loeb and Barkana 2001; Sheth and Tormen 
2002, Scannapieco and Barkana 2002). 

Some statistical characteristics of the LSS has been de- 
rived in our previous papers (Demiariski & Doroshkevich 
1999, hereafter DD99; DDMT). They are based on the 
nonlinear theory of gravitational instability (Zcl'dovich 
1970; Shandarinfe Zel'dovich 1989) applied to a CDM- 
like broad band initial power spectrum of perturbations. 
This approach allows to outline the general tendencies of 
the LSS evolution and demonstrates a leading role of the 
initial velocity field and of the successive merging of struc- 
ture elements in formation of the observed and simulated 
large scale matter distribution. It links some quantitative 
characteristics of galaxy walls with the initial power spec- 
trum. 

Comparison of the theoretical expectations with mea- 
sured statistical parameters of the most conspicuous wall- 
like component of structure was performed in DDMT 
and in Doroshkevich et al. (2002) for the mock 2dF sur- 
vey (Cole et al. 1998), the SDSS Early Data Release 
(Stoughton et al. 2002), the Las Campanas Rcdshift 
Survey (Shectman et al. 1996) and the Durham/UKST 
Redshift Survey (Ratcliffc et al. 1998). This analysis con- 
firms that the walls are gravitationally confined and partly 
relaxed Zel'dovich pancakes formed presumably due to the 
ID collapse of matter. Such interpretation was already 
proposed in Thompson & Gregory (1978) and Oort (1983) 
just after observations of the first walls. 

From this comparison it follows that the main mea- 
sured characteristics of walls are quite well expressed in 
terms of the time scale, the coherent length and correla- 
tion functions of the initial velocity field, set by the power 
spectrum of initial perturbations. The time scale and the 
coherent length are expressed through the spectral mo- 
ments and amplitude of initial perturbations, and the ba- 
sic parameters of the cosmological model. This analysis 
provides independent estimates of the amplitude of ini- 
tial perturbations with a scatter <20% and verifies their 
Gaussianity. 

In this paper we extend the model considered in DD99 
and DDMT by taking into account the action of factors 
responsible for the evolution of structure elements after 
their formation. Among the most important are: the ex- 
pansion and compression of LSS elements in the transverse 
directions, the small scale damping of initial perturbations 
caused by the random motions of DM particles and Jeans 
damping, and the acceleration of cloud formation within 
larger structure elements, what creates a large scale bias 
between the spatial distribution of DM and luminous mat- 
ter. As before, this analysis is based on the Zel'dovich the- 
ory of gravitational instability which describes the process 



of LSS formation as a successive compression of matter 
along three orthogonal directions. This implies a succes- 
sive transformation of pancakes into filaments and fila- 
ments into clouds with a progressive growth of the masses 
and sizes of structure elements due to their merging. 

First of all, this approach allows one to estimate ana- 
lytically the rates of formation of LSS elements and their 
mass functions in a wide range of redshifts. For all LSS ele- 
ments - pancakes, filaments and high density clouds - the 
mass functions are found to be quite similar to each other. 
They are also similar to those predicted by the Press- 
Schechter formalism and its extension for elliptical clouds 
(Sheth and Thormen 2002). This similarity demonstrates 
the generic dependence of the characteristics of the LSS 
on the initial power spectrum and indicates that the shape 
of collapsed clouds influences the rate of collapse but does 
not change significantly the mass functions. These results 
complement the Press-Schechter approach and allow one 
to obtain independent estimates of the mass functions for 
a wider class of objects. 

A special problem is the interaction of small and large 
scale perturbations. As is well known, it accelerates or 
decelerates formation of clouds within compressed or ex- 
panded LSS elements. In simulations this interaction man- 
ifests itself by a strong concentration of halos within fila- 
ments and walls at all redshifts. In observations it is seen 
as a strong concentration of galaxies within the richer fila- 
ments and walls. In contrast, a significant fraction of both 
DM and baryons remains within low mass pancakes which 
are more homogeneously distributed than the luminous 
matter. The approach used in this paper allows one to 
quantify this interaction. We show that the formation of 
high density halos and galaxies is modulated by the large 
scale initial velocity field what explains qualitatively the 
large scale bias. 

Using our approach we can describe the basic proper- 
ties of long-lived pancakes observed as the Ly-cv forest in 
a wide range of redshifts. In this paper we discuss the evo- 
lution of two most important characteristics of such pan- 
cakes: their mass function and their mean number density 
along the line of sight. Both characteristics are evidently 
changing with redshift due to the formation of new and 
merging of old pancakes, and their transverse compression 
and/or expansion. These characteristics depend also upon 
the cutoff of the initial power spectrum at small scale. 

Both characteristics play a key role in the discussion 
of observed evolution of the Ly-cv forest and we used 
them (Demiariski, Doroshkevich & Turchaninov 2003) to 
show that the observed properties of absorbers are consis- 
tent with theoretical expectations. These results restrict 
also the spectral moments and the mass of dark matter 
(DM) particles to Mdm > 1-5 keV. This approach can be 
also used to measure the shape of initial power spectrum 
down to scales of <~ 30/i _1 kpc (Demiariski & Doroshkevich 
2003). 

This paper is organized as follows: In Sees. 2 basic re- 
lations are introduced. In Sec. 3 the redshift dependence 
of the expected matter fraction assigned to various types 
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of structure elements is found. In Sec. 4 we discuss the in- 
teraction of large and small scale perturbations. In Sec. 5 
the joint mass functions of DM structure elements are con- 
sidered. In Sees. 6 & 7 we discuss the redshift evolution of 
statistical characteristics of filaments and pancakes. Short 
conclusion can be found in Sec. 8. Some technical details 
are given in Appendixes A & B. 

2. Basic statistical characteristics of the 
Zel'dovich theory 

In this section we present the basic statistical charac- 
teristics of Zel'dovich approximate nonlinear theory of 
gravitational instability used later to describe the pro- 
cess of structure formation and evolution. Main ideas and 
characteristics were already introduced in DD99 and are 
repeated here without discussion. Some definitions are 
improved and corrected what makes the approach more 
transparent. 

In the Zel'dovich theory (Zel'dovich 1970; Shandarin 
& Zel'dovich 1989) the Eulerian, rj, and the Lagrangian, 
qi, coordinates of particles (fluid elements) are related by 



n = (1 + z)- x [qi - B(2)S<(q)] 



(1) 



where z denotes the redshift, B{z) describes growth of 
perturbations in the linear theory, and the gradient vec- 
tor Sj(q) = d<fi/dq~i characterizes the spatial distribution 
of perturbations. The Lagrangian coordinates of a parti- 
cle, qi, are its unperturbed coordinates in the real space, 
ri(z = 0) = qi. For the spatially flat ACDM cosmological 
model the function B(z) can be approximated (DD99) as 
follows: 

fl A + 2.2fi ro (l + zf 



VIa + TO = 1 



(2) 



1 + 1.2fi m 

Here Q m &zQ\ arc the dimensionless matter density and 
cosmological constant. This fit is reasonably accurate for 
all z and is normalized at z = by the condition B(0) = 1. 

In this paper we consider only power spectra with the 
Harrison - Zel'dovich asymptotic, p(k) oc k, at k — > 0, 
and CDM-like or WDM-like transfer functions, T 2 (fc), in- 
troduced in Bardeen et al. (1986, hereafter BBKS): 



p(k) 



A 2 k 

47Tfc4 



f- ) , k = a^Mpc" 1 
ko, 



(3) 



where A and h = -Ho/100 km/s/Mpc are the dimensionless 
amplitude of perturbations and the Hubble constant, k is 
the comoving wave number. The same approach can be 
applied for arbitrary initial power spectra. 

2.1. Coherent lengths and correlation functions of 
initial density and velocity fields 

For the spectra @ , the coherent lengths of initial density 
and velocity (or displacement) fields, l p and l v , are ex- 
pressed through the spectral moments, namely, to_2 and 
m , (DD99): 



m 



x 3 T 2 (x)D w (x)Dj(x)dx, 



m_ 2 



xT 2 (x)dx fa 0.023, x = k/k . 



fc 0% /m_ 2 



6.6 



Mpc, l p = q l v = 



r 3/2 

ko m 



(4) 



Here Dw and Dj describe damping of perturbations 
caused by the random motions of DM particles (see 
BBKS) and suppression of formation of baryonic clouds 
due to gaseous pressure (Jeans damping, see, e.g., 
Matarrese & Mohayaee 2002), respectively. Analysis of 
simulations and observed Ly-a forest (Narayanan et al. 
2000; Barkana, Haiman & Ostriker 2001; Demianski, 
Doroshkevich & Turchaninov 2003) allows to estimate qo 
and Mum as follows: 

Mdm > 1 - 5keV, q ~ (0.5 - 1) ■ 10~ 2 . (5) 

The characteristic masses of DM clouds associated with 
the coherent lengths l v and l p are 



M„ 



An 



14 



<p>K 



2 ■ 10 



: M Q , M p = qlM v . 



(6) 



Large difference between the coherent lengths l v & l p 
indicates that formation of galaxy filaments and walls ob- 
served at small redshifts is mainly driven by the larger 
coherent length, l v . The small scale correlations of density 
field characterized by the coherent length l p are more im- 
portant during formation, at high redshifts, of low mass 
objects such as the first population of galaxies and Ly-a 
absorbers. 

For spectra under consideration, the normalized corre- 
lation functions of displacement and density fields, 

A 2 m_ 2 , , 
(7) 



o ((q-S(q 1 ))(q.S(q 2 ))) 2 
€v{q) =3 — , cr 



k 2 



o p = A mo 



q=(qi-q2)Au, g=|qi-q2|/i«, 

were approximated in DD99. Here, as before, qi & q2 are 
real unperturbed coordinates of two particles at z — 0, tr 2 
and a 2 are variances of the displacement and the density 
fields, and A is the amplitude of initial perturbations in- 
troduced in © . Due to homogeneity and isotropy of both 
the background matter distribution and perturbations the 
correlation functions depend only on the differences of co- 
ordinates |q|. For the most interesting case q <C 1, q < 1, 
the correlation functions can be written as follows: 

* ~ 1 g 2 qo + Zy/q 2 , * 

Vq 2 + qo 9o + vV + 
f _qo_ id£v ^ qo 

15 q*dq q dq ~ 



3 

ql 



q 



for q < qo < 1 , (9) 



and ^ < 1, for qo < q. At q -C qo both functions, £ v 8z£ p , 
only weakly depend upon the higher spectral moments, 
77i2, TO4, .... Some of the more cumbersome correlation and 
structure functions of perturbations are discussed in DD99 
and are presented in the Appendix A. 
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2.2. Basic relations 

As was shown in DD99, to describe the formation of struc- 
ture elements the basic equation QJ has to be rewritten 
using the differences of particle coordinates and displace- 
ments. The separation of two particles with Lagrangian 
and Euler coordinates qi&q2, and ri&r 2 , respectively, 
is described by the equations (i=l,2,3): 



An = (n - r 2 )j = 



I, 



1 



111 



(10) 



ASi = Si(qi/I„) - Si(fln/l v ), q % = 



(qi ~ qa)i 



Evidently, for Gaussian initial perturbations, the PDF 
of ASi is also Gaussian with 



(A5,fe)>-0, -^(AS*(qi)) 



2[1 (11) 



These relations show the symmetry between expansion 
and compression of matter and the progressive growth of 
(ASffe)) oc \q t \ for larger (1 > \q t \ > q ). These ex- 
pressions indicate also a relatively slower rate of evolution 
of objects with larger |^|, because it is defined by the ratio 
ASi/\qi\ oc a q /q cx q^ 1 / 2 . 

According to the relations l|10(l . when two particles 
with different Lagrangian coordinates qi and q 2 meet at 
the same Eulerian point r a caustic - Zel'dovich pancake 
- with the surface mass density < p > |qi — q 2 | forms. 
Following DD99 we assume that all particles situated be- 
tween these two boundary particles are also incorporated 
into the same pancake. This assumption is also used in 
the adhesion approach (see, e.g., Shandarin & Zel'dovich 
1989). Comparison of statistical characteristics of pan- 
cakes with simulations (DD99; DDMT) verifies this ap- 
proach and shows that long lived and partly relaxed richer 
walls accumulate a significant fraction of compressed mat- 
ter. The same condition can be used to describe formation 
of filaments and clouds in the course of successive collapse 
along two and three axes. 

For the general description of structure formation it is 
convenient to combine the three equations (|10[) into one 
scalar equation, namely, 



q%An 



[l-r(z)F(q)], 



(12) 



T(z) = ^-B( Z )=r B( Z ), nq) = ^|M. (13) 
V3l v T s q z 

The random function F(q) with the Gaussian PDF and 

(F)=0, (F 2 )=<r|(q), (14) 

characterizes the evolution of a region defined by the given 
limits of q. The dispersion a p (q) is expressed through the 
structure functions introduced in Appendix A and DD99. 
The conditions F(q) > 1 and F(q) < separate the col- 
lapsed and expanded parts of the volume, the condition 



1 > F(q) > separates part of the volume that will col- 
lapse later. This approach allows to obtain main charac- 
teristics of structure. However, its practical applications 
are not simple due to complicated form of the function 
<J F (q). 

2.3. Characteristics of the deformation field 

In DD99 we assumed that the deformation field is domi- 
nated by two lowest harmonics what allows to character- 
ize the formation of structure by three weakly correlated 
components of the displacements. Analysis performed in 
Appendix B shows that, for the spectra J3J), this assump- 
tion is valid with a precision better than 10%. For the most 
interesting cases qo <C q <Cl and q <C qo *Cl, the normal- 
ized amplitudes of several lowest spherical harmonics of 
the deformation field, bf, are given by: 



« 0.533, 



0.55, 



b\ 



0.381, 



bl « 0.44 > bl, 



0.037, 



bl. 



0.014,. 



and the contribution of higher order harmonics with I > 
4 is only ~ 1% for low mass clouds with q < q <C 1, and 
it reaches ~ 8% for more massive clouds with go *C q < 1. 

These results justify the assumption made in DD99 
to neglect higher order harmonics of perturbations with 
I > 4 and they confirm that the formation of struc- 
ture can be approximately described by the spherical and 
quadrupole components of the deformation field. The in- 
fluence of higher harmonics, even with small amplitude, 
leads to a small scale disruption of the compressed clouds 
because of the strong instability of thin pancake-like con- 
densations and, so, to the formation of internal structure 
of clouds. 

These results indicate also that when the process of 
structure formation is described by the function F(q) the 
ellipsoidal or, in the case of two dimensional problem, el- 
liptical volumes/areas are preferable. 

2.4. Simple approximation 

When only two spherical harmonics are taken into ac- 
count, the general deformation of any cloud can be de- 
scribed by its deformations along the three orthogonal 
principal axes, namely, x\, x<i, &X3. In this case, we can use 
a simpler approach and consider again the three equations 
(|10|l instead of 1)12(1. Applying this approach it is possible 
to obtain some approximate characteristics of the LSS. 
However, its abilities are restricted and some important 
problems can be solved only with the general approach 

m- 

The numbering of principal axes is arbitrary but, fur- 
ther on, we will usually assume that 



ASi/gi > AS 2 /q 2 > AS 3 /q 3 



(15) 



and the cloud collapses fastest along the first axis whereas 
slower collapse - or even expansion takes place along the 
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third axis. For q t — ► this choice agrees with the order- 
ing of principal axes of the deformation tensor (Zel'dovich 
1970; Shandarin & Zel'dovich 1989). 

The correlations of differences of displacements along 
the principal axes are relatively small (DD99), 



3(AS i ( gi )AS J -(g J -)) < 1 



(16) 



what allows us to consider, in many problems, the compo- 
nents of the deformation field as uncorrelated. Dispersions 
a s & a q were defined in (|7lfc HT|l . 

For quantitative estimates it is convenient to rewrite 
1101) in a dimensionless form 



Ar, : 



1 + z 



r(z)si 
A*(<7i) 



Si 



(17) 



<j q y/2$ - e„fe)] 



, 77(^,2:) = —=—-. (18) 



\/2r(z) 



The probability distribution function (PDF) of dimension- 
less differences of displacements Si is Gaussian and it is 
given by: 

d T (s ;'tr S3) = «p[-o-5(-? + ^ + . (19) 

Here the factor 6 appears because of the imposed restric- 
tions (f 1 51) on the range of variables. 

As is seen from (|S1 IT%|) . for two limiting cases (ft » go 
and qi <C go, we have 



K?)«-2~. ^=27V2' 



(20) 



T> V = Vo = 9<<7o<l, 21 

3 tV 6 



and the two particles under consideration moving along 
the i th axis will cross when s, > V^r/i- For smaller (ft < 
go 1 the condition of the crossing does not depend on 
the size of the pancake, it means that collapse of all small 
pancakes takes the same amount of time. 

The parameter go 13' characterizes the damping 
scale in the initial power spectrum and discriminates re- 
gions of strong and small correlations of the initial den- 
sity and velocity fields. As is seen from (|21() . the condition 
7/0 = 1 introduces also the typical redshift, z — z r , 



{ \ / *° 

T(Z r ) =T r = 4/—, 



1+2 



1 + 1.2Q, 

2 . 2 ^1 ™ 



1/3 



(22) 



At larger redshifts, z > z r , r < r r , the formation and evo- 
lution of the LSS elements strongly depends on the small 
scale correlations, while later on their influence progres- 
sively decrease. For fi m = 0.3, go = 10~ 2 and t ~ 0.3 we 
have z r « 6.5 what coinsides with the observed low limit 
of the redshift of reionization of the Universe (Fan et al. 
2001). 



3. Fraction of matter accumulated by structure 
elements 

The fractions of matter accumulated by pancakes and fila- 
ments were found in DD99. Here we extend this approach 
and introduce three cumulative distribution functions for 
more detailed description of fractions of matter accumu- 
lated by pancakes, filaments and walls with different sizes 
and rate of expansion and/or compression. In turn, this 
extension allows one to obtain the mass functions of all 
three kinds of LSS elements and to take into account the 
further evolution of pancakes after their formation. These 
problems are discussed in Sees. 5 & 7. 

3.1. Cumulative distribution functions of displacement 
differences 

As is seen from (|18f) . the conditions Ari < and si/^/2 > 
V^Qi) — Vi are equivalent and both define collapse of a 
cloud along the axis of the most rapid compression. To dis- 
criminate between pancakes, filaments and walls, we have 
to use also conditions which restrict the motion along two 
other directions. For more detailed description of pancakes 
we have also to restrict the rate of pancake expansion. 
These conditions are consistent with restrictions ()15|) and 
they can be summarized as follows: 

^>C>, < 3 >^>-,(23) 

where the parameters Q restrict the rate of deformation 
in different directions. 

The first condition in (|23|l assumes that for a given r(z) 
the clouds with £1 = 771 > 770, n(qi) > %/2~Ci T are collapsed 
along the first axis. For ( 2 = ^(92) > Vo, C3 = r)(q 3 ) > 770, 
the next two conditions exclude from the consideration 
collapsed filaments with /i(c?2) > V2C2T and clouds with 
A* (93) > \/2C3 r - F° r C3 < C2 < they restrict the expan- 
sion rate of expanding pancakes. 

Similar conditions, namely, 

(25) 

with ^1 = 771 > C2 = V2 > ^o, restrict the minimal sizes of 
filaments and, for Ci = Vi > C2 = m > C3 = V3 > Vo = 
n(q = 0), of clouds collapsed at a given t(z) . 

Integration of the PDF l|19|) under conditions ll'-'il 
I24I &I25 J1 results in the cumulative PDFs characterizing the 
evolution of the LSS elements: 

W P (Ci, C2, Ca) = g(l ~ er)(l + e 3 ))(l + 2e 2 - e 3 ) , (26) 
W>(Ci, Ca, Cs) = |(1 - ei)(l + e 3 )(l + e x - 2e 2 ) , (27) 



W c /(Ci,C2,C3) = -(l-ei)[l + ei + e ;-3e 3 + 



(28) 
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3e 3 (e 3 - ei) - 3(e 2 - e 3 ) 2 ] , 

where ei=erf(d). 

Functions W p , Wf & W c i describe the fractions of mat- 
ter assigned to pancakes, filaments and clouds, respec- 
tively, under conditions (gH For C2 = C3 = 
i]2, e 2 = e 3 , expressions (|21fllT7|l become identical to those 
obtained in DD99. 

3.2. Maximal fraction of matter accumulated by 
structure elements 

From the relations 1271 it follows that in the 
Zel'dovich theory the maximal fractions of matter which 
can be accumulated by clouds, filaments and pancakes, for 
q t = 0, t]i = rjo = r^ 1 (z)^q /6, are 

Woi « ierfc 3 (?7 ), W f = ^crfc 2 (r, )(l + erf (%)) , (29) 
V^^erfc( % )(l+erf(7/o)) 2 . 

We see that during the early evolutionary stages, when 
z ^ z r , t <C r r , rjo 3> 1, the fractions 1)29(1 increase 
cx exp(— ki]q) with k — 1,2,3 for pancakes, filaments 
and clouds. During later evolutionary stages, when z <C 
z r , t ^> r r , 770 <C 1, we get from (I2H1) : 

Wd « 0.125, W f kW p k 0.375 . 

In the Zel'dovich theory motion of matter along the 
three principal directions as described by the equation 
<|10l) is considered to be independent. It means that the 
interaction of matter is neglected and therefore in this 
approach the matter fraction accumulated by filaments 
and clouds and the rate of collapse are usually underesti- 
mated. For small r, the acceleration of the collapse along 
the xi and x 3 axes after the collapse along the x\ axis 
and the interaction of small and large scale perturbations 
increase the fraction of matter accumulated by the col- 
lapsing cloud by a factor of about 2. The same effects 
also shift the boundary between structure elements and 
increase the final fraction of matter assigned to filaments 
and clouds. The popular Press-Schechter approach faces 
the same problems because the assumption of spherical 
collapse overestimates the collapse rate and the fraction 
of collapsed matter. 

The cross correlations of orthogonal displacements as 
given by (|lfi|) . though small, increase the maximal frac- 
tion of matter accumulated by clouds and filaments up to 
W c i ~ 0.2 and Wf ~ 0.3, respectively, and decrease the 
fraction of matter accumulated by pancakes to W p ~ 0.3. 
These maximal fractions are plotted in Fig.^versus r/ 2 /r/ 2 
for all three types of structure elements. As it is apparent 
from Fig. Q already at r = r r , r\ = 770, ~ 25%, 8% and 
2% of matter is accumulated by pancakes, filaments and 
clouds, respectively. At t < T r , these mass fractions in- 
crease cx exp(— ray 2 ), k = 1, 2, 3 what coincides with (|29|l . 

These estimates of the mass fractions are quite reason- 
able for large redshifts, when the formation of structure 
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Fig. 1. Fraction of the DM component accumulated by 
pancakes and walls (solid line), filaments (dashed line) 
and clouds (long dashed line) plotted versus rf 2 . These 
fractions at the moment 77 = 1, r = r r , are shown by thin 
straight lines. 

elements from the dispersed matter dominates. At small 
redshifts the merging of the LSS elements distorts these 
estimates and, for example, some fraction of filaments and 
clouds is incorporated into richer walls. These estimates 
can be improved by a comparison with high resolution 
simulations and observations of Ly-a forest at high red- 
shifts. 

3.3. Distribution function of the sizes of structure 
elements 

The relations l|26l 1271 128fl are monotonic functions of Q 
and, for Q — rji > 7/0, they represent the cumulative 
three parameters distribution functions of structure ele- 
ments with respect of their sizes measured by differences 
of Lagrangian coordinates, qi,q2,q 3 , for particles bound- 
ing these LSS elements along the principal directions. This 
joint PDF for all structure elements is 

N{qi,q 2 ,q3) oc exp(-7y 1 - r/ 2 - r} 3 ) — — — (30) 

dqi dq 2 dq 3 

VO < Vi < OO) < 1i < 00 • 

Similar relation allowing for the correlations (|16|) can be 
obtained in the same manner. 

The expression l(30|l does not discriminate filaments 
and pancakes with respect to their expansion rates de- 
scribed in l|2()l) and 127f) by negative £2 and £ 3 . To al- 
low such discrimination it is necessary to link (2 < and 
£ 3 < with the size of structure elements. The symmetry 
of the basic characteristics of structure elements with re- 
spect to positive and negative ASi discussed in Sec. 2.2 
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and the relations (|llf> indicate that such an extension can 
be achieved by imposing the following limits in (|30|) : 



< qi < oo, — oo < q2 < oo, — oo < q% < oo . 

< qi < oo, < qi < oo, — oo < 93 < 00 , 
< qi < 00, < qi < 00, < 93 < 00 , 



(31) 



High resolution simulations show that at high redshifts 
substantial fraction of filaments accumulate high density 
halos. In observed redshift surveys galaxies are mainly sit- 
uated within filaments and walls and the population of 
isolated galaxies is quite small. This strong galaxy con- 
centration within the LSS elements - large scale bias - 
can be naturally explained by the interaction of small and 
large scale perturbations responsible for the formation of 
galaxies, filaments and walls, respectively. 

As was noted in DD99, the statistical approach allows 
to describe this interaction. Here we are returning to this 
problem and propose a simple quantitative measure that 
illustrates effects of this interaction. We compare the pro- 
cess of formation of two identical pancakes with a surface 
density mi — (p)qi, and /Ji = /j(9i). We assume that 
the first - a reference - pancake is formed at the 'time' 
moment ti(z\) and is characterized by the function 



V2- 



(l + *l) 



y/2- 



as given by (|18l I21[) . For comparison, let us consider the 
process of formation of the second pancake with the same 
size, q2 = qt, /i(9 2 ) = ^2 = Mi: at redshift Z2 and the 
successive formation, in the same region, of a pancake with 
m 3 > m 2 , 93 > 92, M93) = > A*2 at redshift z = z 3 < 
Z2, and at the 'time' moment T3{zz) > t 2 (z2). We assume 
that the larger pancake with the size 53 accumulates the 
earlier formed smaller pancake with the size 92 <C 93. 

In this case, the interaction of pancakes can be de- 
scribed by replacing 772 in 1)261 130(1 with the effective pa- 
rameter r] e ff defined by the relation (DD99): 



Veff = {V2 ~ r23?73)/ \Jl ~ ^3, 
2/X2M3 



(32) 



^23 



9293 



92 - 93 



92 + 93 



If both pancakes are small, 1 3> qo > 93 S> 92, then 
^23 — * 1, and the strong small scale correlation of density 
field decreases the survival probability of smaller pancakes 
and both pancakes are actually formed at similar redshifts 
Z2 ~ 23. However, later formation of a larger pancake with 



».o r 



for pancakes, filaments and clouds, respectively. Here neg- 
ative Q2,93 are associated with expanding filaments and q j 
pancakes. q 1 

4. Interaction of small and large scale 
perturbations. Large scale bias 
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Fig. 2. Fractions of matter, f m , accumulated by structure 
elements situated within HDRs (solid lines) and LDRs 
(dashed lines) are plotted vs. the threshold overdensity 
for simulated DM distribution at three redshifts. 

1 ^ 93 > 9o, 93 92, M3 /-*2, accelerates compression 
of the smaller pancake, it increases the effective ampli- 
tude of perturbations and shifts formation of pancakes to 
larger redshifts. For T23 ~ M2/M3 "C 1, and at larger red- 
shifts, when rji oc 1 + z, the general expression 13 2|) can be 
rewritten more transparently as follows: 



Veff « % 1 



l + z 3 



(Z2 ~ 23) ■ 



(33) 



1 + Z 2 ) V^To 

Statistical characteristics of pancakes formed at z = Z\ 
and z = Z2 will be the same for the same value of 



Veff = 771 = (1 + zi) 



M9i) 
V2r 



However, as is seen from H33fl . this occurs already at the 
redshift 



Z2 = Zi 



1 + z 3 > zi + 1 



This relation gives the quantitative estimate of the in- 
fluence of interaction between the large and small scale 
perturbations, which increases the actual redshift of pan- 
cake formation. Pancakes with the same size qi are now 
formed earlier at redshift z% > Zi, instead of at z\. 

The suppression of pancake formation within expand- 
ing regions can be considered in the same manner. Thus, 
for example, for negative ASi and — Si > \/2r](qi,z) we 
have instead of l|33|l 



Veff ~ V2 1 + 



1 



1 



Z2 



M(92 



(2 + z 2 + z 3 ) , 



\/2to 

and from the condition rj e f f = Vi we S e t 
Z3 < zi 



(34) 



Z2 



zi-1 



1 
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what illustrates the influence of this interaction on the 
moment of formation of pancakes. 

The influence of large scale perturbations on the pro- 
cess of formation of small scale objects can be considered 
as the manifestation of large scale bias. Fraction of matter 
accumulated within high density clouds increases rapidly 
with time (see Fig.^) and halos formed at larger z contain 
only a small fraction of mass of halos formed at smaller 
z. But this bias increases the redshift of subclouds forma- 
tion and their densities, promotes the transformation of 
DM clouds into observed galaxies, and makes the internal 
structure of clouds more complex. Regular distribution of 
young galaxies within high density filaments suppresses 
their feedback and fosters further formation of galaxies. 

Direct calculation shows that, for larger clouds with 
the correlation coefficient r 12 <C 1, the influence of the 
nearest collapsed cloud dominates, and the generalization 
of (|521 fe l55|) to the case of multi step collapse at z 2 > z 3 > 
za > ■■■ practically does not change the dominant role of 
the first term. This means that the multi step collapse has 
to be considered step by step taking into account only the 
next generation of high density condensations. This means 
also that such interactions can be described by the theory 
of Markov processes. 

4.1. Comparison with simulations 

To compare these expectations with simulations we 
use the high resolution simulation (Klypin et al. 1999; 
Schmalzing et al. 1999) performed with adaptive code in a 
box of (eO/i^Mpc) 3 with N to t = (256) 3 particles for the 
Harrison-Zel'dovich primordial power spectrum and the 
BBKS transfer function. This simulation approximately 
corresponds to qo ~ 10~ 2 , Mum ~ 1 keV. The matter 
distribution at redshifts z = 4, z = 1 and z = were 
analyzed. 

To test the impact of environment on the properties 
of high density clouds the full samples were separated 
into subsamples of high and low density regions (HDRs 
and LDRs). The subsamples of high density regions ac- 
cumulate ~ 40 - 45% of all DM particles within richer 
clouds selected with the friend-of-friend method for the 
threshold density equal to the mean density. Subsamples 
of LDRs were prepared by removing the HDRs from the 
full samples. 

The fraction of compressed matter is plotted in Fig. 
[21 versus the threshold overdensity used for selection of 
structure elements. As is seen from this Figure, at all red- 
shifts the majority of highly compressed matter is situated 
within the HDRs. This excess is moderate at z — and 
progressively increases with redshift. Special test shows 
that the high density clouds selected at z — 4 are equally 
distributed between HDRs and LDRs selected at z = 0. 
This fact is consistent with almost equal concentration 
of galaxies within the HDRs and LDRs at z = in the 
observed surveys. It indicates also that the efficiency of 



the interaction of small and large scale perturbations de- 
creases with redshift. 

5. Mass function and rotation of structure 
elements 

The popular Press-Schechter formalism (see, e.g., review 
in Loeb and Barkana, 2001) focuses main attention on 
the spherical collapse and the achieved final critical over- 
density. It can be extended for elliptical clouds (see, e.g., 
Sheth and Tormen 2002). In spite of so strong restric- 
tions on the process of collapse it successfully describes 
the simulated mass distribution of high density halos. In 
Zel'dovich theory of gravitational instability similar ap- 
proach allows to find the mass function of all high density 
structure elements - clouds, filaments, and walls or pan- 
cakes - without any assumptions about their initial and 
final shapes, stability, relaxation and achieved overdensity. 
In this regard, our approach is much more general then the 
Press-Schechter formalism. 

The close connection of the observed objects with the 
initial power spectrum and the generic origin of galax- 
ies and other observed elements of the LSS is clearly seen 
from the expression (|30|) which applies to all structure ele- 
ments. Due to symmetry of the relation 1|30|1 with respect 
to variables q\ , q 2 & q 3 and because the mass functions 
are normalized, the restrictions (|23l I24lfe l25fl , used for the 
determination of the mass fractions of different structure 
elements, are now of no importance and they lead only 
to renumeration of coordinates. Thus, in the Zel'dovich 
approach, the mass functions of the various LSS elements 
differ only by their corresponding survival probability (see, 
e.g., Peacock & Heavens 1990). 

5.1. Initial shape of collapsed clouds 

To find the mass function of structure elements we rewrite 
the general relation (|30|1 in spherical coordinates 

qi = Rcos8, q 2 — Rsin9 cos (f>, q 3 = R sin 9 sin (j) , 

where < R < 00 and < 9, <\> < ir/2 characterize the 
mass and shape of the collapsed structure elements in the 
Lagrangian coordinates: 

M cx R 3 , tg<f> = q 3 /q 2 , tg9 = sjqj+ql/qi . (35) 

For spherical clouds tg</> = 1, 4> = 7r/4, tg6> = \/2. As 
was noted above, here we will neglect the restrictions Q23I 
F2~Hfc l25|) and differences in limits of qt (|3T|) . 

This result demonstrates that in the Zel'dovich ap- 
proximation the distribution of shapes of initial clouds 
described by the parameters 9 & <fr in (|35|l is continu- 
ous and the finally achieved properties of structure ele- 
ments, such as their shape, energy and overdensity, de- 
pend mainly upon the velocity or deformation field within 
clouds. In particular, even high density compact relaxed 
halos can be formed through the collapse of initially asym- 
metrical clouds. Some information about the initial shape 
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of collapsed clouds is retained in the angular momentum 
of observed clouds discussed below. 

On the other hand, for all clouds, the first step of col- 
lapse is the formation of a pancake-like objects (Zel'dovich 
1970; Shandarin et al. 1995) which are unstable and will 
rapidly break up into a system of low mass subclouds 
(Doroshkevich 1980; Vishniac 1983). This instability stim- 
ulates formation of numerous satellites of the largest cen- 
tral object and makes it difficult to observationally distin- 
guish between isolated galaxies and these satellites. 

For the more massive clouds with 1 > qi ^> go, r)i oc 
y^7, we have from l|30[) : 



d 3 W = 



exp(-i?V(6 l , <^)/8t 2 ) dR 3 dcj)d8 



y/ R 3 cos sin 2<fi 



ip(6, <j>) = cos6 + V2 



sin W cos 



67T 3 / 2 T 3 ' 

- tt/4) . 



(36) 



This relation shows that the fraction of initially massive 
spherical clouds with ip = \/3 is exponentially suppressed 
and asymmetric massive clouds dominate. 

For low mass clouds with qi < go, we get for the mass 
function 



(37) 




d 3 W oc R 3 dR 3 T~ 3 sin 20(1 - cos 26) sin 2<j)ddd<j) , 



and the formation of low mass clouds with M < M p © is 
suppressed. However, in this case the existence of approx- 
imately spherical initial clouds is much more probable. 



5.2. Survival probability of objects 

As it was shown above formation of low mass structure 
elements is usually suppressed because they could be ab- 
sorbed by larger objects formed at the same time and 
in the same region. This process is described by the sur- 
vival probability of structure elements and it is different 
for clouds, filaments and walls. As was shown in DD99, 
for walls formed through the ID compression along the q\ 
axis the survival probability can be taken as 



W, 



rf(vW87 



(38) 



For filaments and clouds formed through the 2D and 3D 
compression the survival probabilities can be taken as 



W, 



^i{^q~^)^i{^fq~^) , 



(39) 



er 



f ( y/ qi /8r 2 ) erf ( V 92 /8r 2 )erf ( ^/q 3 /8r 2 ) , (40) 



respectively. Evidently, so defined survival probability is 
small for low mass objects and W surv — > 1 for massive 
objects. 

Integration of the PDF l|3t)l) corrected for the survival 
probability (|3"51 l3~9l fc ED} . W surv d 3 W , over angular vari- 
ables for given R, go & r allows one to find the mass func- 
tions of structure elements: clouds, filaments, walls and 
pancakes. 




0.01 



(M/<M> 



.1/3 



Fig. 3. Functions N w (thick solid line), Nf (thick dashed 
line) and N c (thick long dashed line) for two values of qo 
plotted vs. (M / {M)) 1 / 3 for r = 0.3. For r « r r = ^fq^J§ 
and r 0.3t>, functions N w are plotted by thin solid and 
dashed lines, respectively. 



5.3. Mass function of structure elements 

The mass function N m and the function M / (M)N m are 
plotted in Figs. H&H for q = 10~ 2 , & q Q = 10~ 3 and for 
the three more interesting values of r = 0. 3, t = r r , & r« 
0.3t>. If t ~ 0.3 corresponds approximately to the present 
epoch (see Sec. 7.1), then r ~ r r (see (j22J) an d r ~ 0.3r r 
describe the early period of structure formation, when the 
influence of small scale correlations of initial density and 
velocity fields is more important. 

As is seen from Fig.|21at r < r r the mass functions for 
all three kinds of structure elements are identical, what is 
a direct consequence of strong correlations of small scale 
perturbations. However, the difference between these mass 
functions increases with time and at r w 0.3 it becomes 
significant, especially for go = 10 -3 . 

As is seen from comparison of Figs. 01 & El in all the 
cases the numerous low mass structure elements contain 
only negligible fraction of mass and for all r the main mass 
is concentrated within structure elements with M ~ (0.2— 
0.7) (M) . For t < r r the impact of the parameter go on the 
shape of mass functions is negligible and, for both chosen 
values of go, the functions MN m (M) are well fitted by 
the Gauss function (|41|) . This result reflects again a strong 
small scale correlation of the initial density field and, as 
was discussed in Sec. 5.2, a small survival probability of 
low mass structure elements. 

At t > T r the mass function becomes wider because 
the continued formation of low mass structure elements 
is accompanied by progressive mass concentration within 
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Fig. 4. Top panel: Functions M/(M) N w (thick solid line) 
and M/{M) N c (thin solid line) plotted vs. (M / (M)) 1 / 3 
for t = 0.3. Fit (|43|1 is drawn by dashed line. Bottom 
panel: functions M/(M) N w for r re t t = y^qo/6 (thick 
solid line) and r re 0.3r r (thin solid line) plotted vs. 
(M/(M)) 1 ^ 3 . Fits pTlfcl4l)l are drawn by thin dashed and 
long dashed lines, respectively. 



V 




10 



0.01 



0.1 



Fig. 5. Functions (M)/t 4 - 3 (triangles, qo — 10~ 2 ) and 
(M)/t 5 - 2 (up down triangles, qo = 10~ 3 ) plotted vs. r. 
Fits l|46l fc l47JI are drawn by solid lines, respectively. 

5.4. Mean mass of objects 

The mean mass of objects is quite sensitive to the coherent 
length of initial density field. For the same qo as above, the 
mean mass of clouds is plotted in Fig. [S] versus r together 
with fits 



o . Tn2 T 4.3 

(M) re _ 9 ^ A M V , for q = 10" 2 



(1 



(46) 



massive elements with M - (M) > M p ©. As is seen m 
Fig. |3 during this period the small scale cutoff of initial 
power spectrum and the parameter qo provide the low 
mass cutoff of the mass function. For such redshifts, the 
mass functions for objects with M < (M) and M > (M) 
are described respectively by power and exponential laws. 

Disregarding the low mass 'tails' of mass functions, we 
fit the mass function for both values of qo by: 

0.4exp(-(x- 0.95)70.17), 



x 3 N m (x) 
x 3 N m (x) 



20x 2 exp(-4.1x)erf(x 2 



0.3r r , (41) 
(42) 



x 3 N m (x) » 12.5a: 2 exp(-3.7a;)erf(x 2 ), r = 0.3, (43) 

where x = (Af/(A/)) 1 / 3 . These fits are plotted in Fig. 03. 
For r re 1, we have 

x 3 N m (x) re S.x 3/2 exp(-3.1a;)erf(a; 9/4 ) . 



This function can be used to characterize the observed 
and simulated filaments and walls which are still forming. 

By definition, the mass functions must satisfy two nor- 
malization conditions 



3 / x N m (x)dx : 
lo 



/•CfO 

3 / x 5 N m (x)dx 
Jo 



1 . 



(45) 



Fits - 0) better satisfy the second condition (gSJ 
while the first one is violated more strongly. These mass 
functions are similar to the expression (|36J) and N rn cx 
exp(— x), for x ^ 1 and r > r r . 



(M) 



3 ■ 10 3 r 5 - 2 

( 1 + T 1.2)5.2 



M„ 



for q = 10 



(47) 



which describe quite well the redshift evolution of (M) for 
r > r r . For t < r r the much slower evolution of the mean 
mass is described by 



(M) re t 3 M v , for q Q = 10" 



(M) re 4 • IQ- 2 t 3 M v , for q Q = 10" 



(48) 



(49) 



what again emphasizes the impact of the small survival 
probability of clouds at r < r r . 

For r > r r , the mean masses of pancakes and filaments 
are less than those for clouds by a factor of ~ 1.5 - 2. 



(44) 5.5. Comparison with the Press - Schechter formalism 



The mass function N m describes all structure elements 
- clouds, filaments, and walls or pancakes - without as- 
sumptions about their shapes and achieved overdensity. 
However, the relation (|44|l is quite similar to the Press- 
Schechter mass function for scale-free power spectra and 
k corresponding to typical objects, p cx fc -2 : 



M 

W) 



Npsdi = 



45Vtt 



e i / 6 ex P (-e i/a ) 



1/3 N 



(50) 



f= 1.875 M/(M), (0 = 15/8. 
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Fig. 6. Mass functions of galaxy clouds, N mem / (N mem ))N m , selected at redshifts z — (left panels) and 
z — 4 (right panels) in HDRs and LDRs for two threshold linking lengths. Fits 143J1 and (|52|l are plotted 
by solid blue and green lines, respectively. 



In turn, for larger redshifts, r ~ r r , and for a power 
spectrum with a cutoff at k ~ fc ma; , both approaches pre- 
dict the suppression of formation of low mass objects with 
M < (M) (Loeb & Barkana 2001). This similarity is quite 
apparent as both relations are based on the same initial 
power spectrum. It indicates that the difference between 
the Press-Schechter and a more general Zel'dovich ap- 
proach is quantitative rather then qualitative and these 
approaches are complementary to each other. 

5.6. Impact of coagulation processes 

All mass functions discussed above are related to the pro- 
cess of formation of structure elements and they do not 
take into account the later nonlinear evolution described 
by the coagulation equation (Smoluchowski 1916; Silk & 
White 1978). This is not so important for walls and fil- 
aments for which the merging and coagulation are con- 
trolled mainly by the initial velocity field. But the non 
linear evolution can essentially distort the mass function 
of high density clouds accumulated by richer walls and 
filaments. 

In contrast with expressions (|4 11441) and (|50() the co- 
agulation process leads to the mass function 



N c (x) oc x~ v exp(-ar), x = M/(M) > x 7J 



(51) 



where the power index v ~ 3/2 depends upon the aggrega- 
tion rate (see, e.g., Silk & White 1978). This mass function 
strongly differs from those discussed above. This means 
that the shape of observed and simulated mass functions 
measures the possible influence of coagulation processes. 



5. 7. Comparison with simulations 

The theoretical fits (|43l & I50fl can be compared with the 
mass functions found for matter distributions at redshifts 
z = 4, and z — in the high resolution simulation (Klypin 
et al. 1999; Schmalzing et al. 1999) discussed in Sec. 4.1. 
As was described in Sec. 4.1, at both redshifts the full sam- 
ples were divided into subsamples of high and low density 
regions (HDRs and LDRs) dominated by walls and fila- 
ments. 

The mass functions of high density clouds found for 
the HDRs and LDRs separately are plotted in Fig. The 
basic parameters of the same samples of selected clouds 
are listed in Table 1, where Sthr and (5) are threshold and 
mean overdensities of clouds above the mean density of 
the sample, f pn t is the fraction of points accumulated by 
clouds, N c i and (N mem ) are the number and the mean 
richness of clouds. At both redshifts, very massive struc- 
ture elements are formed through the percolation process 
and they cannot be described by the expressions H43l fel5T7 | > ■ 
For this reason, they were excluded from the analysis of 
the mass distribution. However, these clouds are included 
in estimates of the matter fraction f pnt accumulated by 
structure elements as they present actual filaments and 
walls. The difference between both f pnt and (N mem ) for 
clouds selected with the same S t hr within HDRs and LDRs 
confirms a significant impact of environment on the prop- 
erties of clouds. 

At both redshifts, the samples selected with high 
Sthr represent properties of small fraction of high density 
clouds, while samples selected with small S t hr are formed 
mainly by unrelaxed filaments and walls. The cutoff of 
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Table 1. Parameters of structure elements selected in 
HDRs and LDRs at z=0 & 4. 





Sthr 


(5) 


fpnt 


N cl 


{-^V mem } 








2=4 






HDR 


490 


3 007 


0.015 


2 187 


101 


HDR 


1.6 


15 


0.37 


5 311 


313 


LDR 


7.1 


139 


0.004 


963 


75 


LDR 


1.6 


15 


0.036 


6 104 


108 






z = 








HDR 


86 


1 200 


0.18 


11 155 


273 


HDR 


22 


645 


0.33 


2 908 


342 


LDR 


22 


712 


0.21 


8 484 


270 


LDR 


1. 


62 


0.43 


9 629 


261 



The few richest clouds formed through percolation process 
are excluded from calculations of (N mem ). 



the simulated mass functions at low masses caused by 
a finite resolution increases the measured mean richness 
of selected clouds as compared with theoretical expecta- 
tions 14611 . Because of this, the mean richness was used 
as a fit parameter in the relations l|43|) plotted in Fig. 
However, the shape of the mass function (|43[) was not al- 
tered. As is seen from Fig. this expression fits well the 
simulated mass functions for both high density clouds and 
only partly relaxed filaments and walls. 

However, at the redshift 2 = the simulated mass 
functions are well fitted also by a power law 

xN m (x) oc a;" ' 85 , x = N mem /(N mem ) , (52) 

similar to (|51|l . It indicates the possible influence of the 
coagulation processes on the parameters of selected high 
density clouds and even filaments and walls. At the red- 
shift z = 4 this influence is weak and the mass functions 
are evidently exponential. 

It is especially important that at both redshifts the 
relations 1(431 fc l5T) | successfully reproduce the mass func- 
tions of unrelaxed filaments and walls which are far from 
the spherical shape. This fact demonstrates a moderate 
influence of the shape of collapsed clouds on their mass 
and the validity of Zel'dovich approach which considers 
the dynamical characteristics of collapsed clouds rather 
than their shape. 



5.8. The angular momentum of collapsed clouds 

Using the statistical approach it is possible to consider 
the angular momentum of the collapsed clouds. For this 
purpose we use the general equation Q together with the 
corresponding expression for the velocity of fluid element 



dr 



H{z)r l -^{P{z)-l)B{z)S l . 
dt 1 + z 



(53) 



where (3{z) = (1 + z)d\nB/dz, H(z) is the Hubble param- 
eter, and B(z) was introduced by As usual, we define 
the angular momentum of a particle as 



ji = djkrjVk — JoeijkdjSk, 



Jo = — (IT? — (54) 



where Cy-fc is the unit antisymmetric tensor and the func- 
tion Jo (z) describes time variations of the angular mo- 
mentum. Let us note that the description of the angu- 
lar momentum through the deformation tensor (see, e.g., 
White 1984) is useful methodically but cannot be applied 
to more massive objects because of the small correlation 
scale of this tensor. 

The angular momentum of a cloud is defined by the 
integral over the corresponding collapsed volume, V, 



(55) 



/(p, q) = 2(pq)G 1 (|p - q|) + (p 2 q 2 - (pq) 2 )G 2 (|p - q|) , 

(the functions G\ & Gi are introduced in DD99 and 
Appendix A). It depends upon the size and shape of the 
collapsed region, statistical characteristics of which are de- 
scribed by l(3u|) . For low mass early formed clouds with 

p<qo,q< qo> we get 



I~ [3(pqr-pV]Azo, 



(56) 



and, for example, for ellipsoidal clouds with axes 
ai, 0,2, & a.3 we have 



/ d 3 pd 3 ql cx — [ofj, + o? 3 + o| 3 ] , 
Jv to 



-a), (57) 



what is identical with the expression found already in 
Doroshkevich (1970). 

For larger clouds with p ^> go, q go, we get 



5(pq) 2 - p 2 q 2 - 2(pq)(p 2 + g 2 ) 
yV + g 2 - 2(pq) 



(58) 



and for spherical clouds (J 2 ) = 0. Numerical integration 
of (55) shows that for large elliptical clouds (J 2 ) depends 
on the axes a; in a similar way as for small clouds 1)57(1. 

6. Statistical characteristics of filaments 

In addition to the fraction of matter accumulated by fil- 
aments (12711 and the mass function of filaments given by 
1(411 I44f) , we will consider here also the linear density of 
galaxies along a filament, T,fu, defined as a mass per unit 
length of filament, and the mean surface density of fil- 
aments, (Jfu, defined as the mean number of filaments 
intersecting a unit area of arbitrary orientation. These 
characteristics depend upon the threshold linking length, 
rink , used for the filament selection, which determines the 
threshold overdensity bounding the filaments. Now only 
richer filaments can be selected in both observed and sim- 
ulated catalogues what restricts its quantitative charac- 
teristics. Because of this, here we will only discuss char- 
acteristics of richer filaments. 

The distribution function of filaments linear density 
describes their frequency distribution with respect to the 
amount of matter per unit length of filaments. As was 
discussed in DD99, this function can be obtained by in- 
tegrating ptljl over all g3 and over the ratio gi/g2- For 
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Fig. 7. Distribution function, Nfu , for the linear den- 
sity of DM particles in filaments selected at three linking 
lengths, rink- Fits (|60|l are plotted by solid lines. 

qo «C 1 , go <C q < 1 and using the survival probability for 
low mass filaments given by (|39[1 we have for the PDF of 
richer filaments with S/y oc g 2 + q\ > < S fu > 

N fa « 7 ^- T exp(- v /3E /i ,/(E /i ,)), (£/«) « 48r 4 .(59) 

In Fig. [7] the PDFs, Afyy, are plotted for filaments 
selected from two simulated DM distributions at z = 
(Jenkins et al. 1998, top panel, Cole et al. 1998, middle 
and bottom panels) for three linking lengths ri nk . The lin- 
ear density of matter, £/,•;, was measured by the ratio of 
the number of points and the length of the minimal span- 
ning tree for each filament of the sample. 

The PDF, Nfa, is well fitted by 



Nfu = a erf 4 [ai(x - x )] exp(-^a 2 (a; - x )) . 
where x = £/«/(£/«) and 



(60) 



a 


= 50, 


«i 


= 2.5, 


«2 


= 27, 


x 


= 0.3, 


o n 


= 180, 


«i 


= 2, 


«2 


= 42, 


x 


= 0.35, 


a 


= 60, 


a i 


= 2.5, 


«2 


= 31, 


xq 


= 0.35, 



for these samples, respectively. Here the cut-off of the 
PDF at x = xq reflects the limited resolution with which 
the filaments were selected. This cut-off increases 
and changes the parameters clq, a%, ai in 1)60(1 with respect 
to the expected ones (|59|) . The measured mean linear den- 
sity, (Efn), depends upon the linking length used for the 
selection of filaments, and for the samples used we have 

rinkTifa = 3.4, rinkT,fu — 3.1, n n feS/« = 3.1 , 



respectively. At the same time, for all samples at 
(£/,:;) the exponential term dominates what is consistent 
with theoretical expectations (|59|l . 

For the surface density of filaments, {(Jfu), and for 



their mean separation, (Dfu) 
tively: 



-1/2 



we get, respec- 



Wf(ri nk )(n) 

(i + ^) 2 (s /i; > 



(Dfi 



W f {r lnk )(n) 



(61) 



Here Wf(ri nk ) and (n) are the fraction of objects accu- 
mulated by filaments selected with a given ri, lk and the 
mean density of objects in the sample. Evidently, Wf <C 1 
for rznfc(n) 1 / 3 < 1 and Wf — ► 0.3 — 0.5 for larger n n k- The 
variations of Wf strongly influence (c/y) and (Dfu) and 
determine how they vary with r;„fc, which is strongly con- 
nected with the overdensity of selected filaments. However, 
these estimates of {erf a) do not take into account possi- 
ble accumulation of filaments by walls and therefore they 
should be corrected by comparing them with observations 
and simulations. 

7. Characteristics of walls and pancakes 

From the PDF of the differences of displacements l|19JI 
it is possible to extract important approximate charac- 
teristics of walls and less massive pancakes which can 
be directly compared with available observations. Some 
of them were introduced in DD99 and successfully com- 
pared with simulated and observed characteristics of walls. 
Other characteristics discussed in this Section are suc- 
cessfully used for description of low mass pancakes ob- 
served as Ly-a absorbers at large redshifts (Demianski, 
Doroshkevich & Turchaninov 2003) and for the deter- 
mination of the initial power spectrum at small scale 
(Demianski and Doroshkevich 2003). 

7.1. Distribution function of wall sizes 

The distribution function of sizes of walls describes their 
frequency distribution with respect to their Lagrangian 
size what is identical to their surface density, m w , defined 
as the mass per unit surface area of the wall at the moment 
of its formation. This distribution function is obtained by 
integrating H30fl over all qi and 03 for go *C 1, q < 1 and 
it is given by: 



1 ®w{qo/q w ', 
|qi - 



2ttt 



■ exp 




(62) 



l v (p) l v 
q w N w (q n 



N w (q w ) dq 



dq w « 8(0.5 + 1/tt)t 2 « 6.55t 2 



®w{y) = ^Vqwdfi(q w ) I dq v 



The factors eri(yJq w /8T 2 ) and 6™ in (|rJ2|l describe the sur- 
vival probability of walls and the influence of the coherent 
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length of initial density field, 90 For q w ^> qo, Q w = l and 
N„, becomes identical to the one found in DD99. 



7.2. Transverse characteristics of walls and pancakes 

For some applications we have to estimate the transverse 
characteristics of pancakes such as the distribution func- 
tion of Lagrangian size and the mean real size of pancakes. 
These characteristics can be found with the method used 
above. 

For the frequency distribution of the walls with 
Lagrangian transverse sizes, 92 & 93 , we get from Ij2()l 1211 



op 

N(rj2)di] 2 



■ exp(-7j|)d7? 2 , Vo < m < 00 . 



(63) 



(r ?2 ) = l/^, (?f 2 ) = (q 2 )/&T 2 = l/2. 



and the same distribution for 773 & (73. In these cases, for 
the low mass objects the merging and percolation are not 
so important and the functions (|63[1 are not corrected for 
the survival probability. 

The mean transverse size of expanded, Ar e , and com- 
pressed, Ar c , pancakes can be found from relations IjlOl 
IT§|) and we have 



(Ar c 



An 



2 I 2tt 



(At 



(At 



7 A 2 

2\ _ '^0 



13AJ 
16 



1+^] , (64) 



10 



13tt 



8t 2 
\ + z 



~ In 



At small redshifts for the ACDM cosmological model with 
T = 0.2, 8tq ps 0.5, the expected mean transverse size of 
walls ~ 0.5^ sa 16ft.~ 1 Mpc is similar to the observed one 
(Doroshkevich et al. 2001). 

7.3. Distribution function of pancakes surface density 

After pancake formation, the transverse compression 
and/or expansion of matter changes its surface density 
and other characteristics. However, the direct analysis in- 
dicates that the PDF of pancakes surface density given by 
H62(l only weakly depends upon these deformations. 

The evolution of the pancake surface area, Span ? is de- 
scribed by the relation p0(l as follows: 



^pan(z) OC 



Ar 2 Ar 3 

9293 



1 



{l + zf 



1 - 



S2 



V2772 



1 - 



S3 



Therefore, for a pancake formed at a redshift z/ with the 
surface density 

Vpan(Zf) = 1 V {1 + Zf) 2 qx , 

the surface density at a redshift z is 

0~pan{z) = &pan{Zf)Sp an (zf)/Sp an (z) . (65) 



This means that the fraction of matter accumulated by 
pancakes with the surface density > a pan depends upon 
their transversal sizes, 92 & 93 , and instead of l|26|) we get 
from 



W a 



erfc 



dx 2 dx3 exp(— x 2 — x\) x 



[1 - x 2 /t] 2 (z)][\ - £3/773(2:)] \ 

[1 - £2/772(2/)] [1 -£3/773(2/)] / 



(66) 



As before, here d z W a jd z qi is the PDF similar to (|30|l and, 
after integration of this PDF over all 92, 93, we obtain l|66|l 
with 92, 93 — > 00. So, we get 



AT^exp^)-^ 

pan 



do, 



■nl = 



'pan 



p 8r 2 (l + z) 2 



If the survival probability of a pancake is erf(77 p ) then 
the normalized PDF of surface density is given by (|62|) . 
where q w and r are replaced by a pan and t z = (1 + z)r, 
and, for example, 



2_2 



(o-pon) =8(0.5 + 1/tt)(1 + zYt 



(67) 



However, if we retain the expression for the survival prob- 
ability used in l|62() then the statistical characteristics of 
pancakes will depend upon both the redshift of pancake's 
formation and its current redshift. 

These results indicate that at high redshifts, when 
(1 + z)t(z) w const, the expansion and compression of 
pancakes approximately compensate each other and the 
PDF (62) or (66), the mean surface density, (a pa n), and 
other average characteristics of pancakes only weakly de- 
pend upon the redshift. They indicate also that, in spite of 
the strong evolution of each individual pancake, the statis- 
tical description 16211 remains valid also when we consider 
each wall and each pancake as formed at their current 
redshifts. 

Application of these results to absorbers observed in 
a wide range of redshifts (Demianski, Doroshkevich & 
Turchaninov 2003) confirms these conclusions. 

7.4. Mean comoving linear number density of walls 

Using the relations H30lfclrl2^l it is also possible to obtain an 
approximate estimate of the mean comoving linear number 
density of recently formed walls, that is the mean number 
of walls per unit distance along a straight line. For richer 
walls with a threshold surface density qthr 3> 9o the small 
scale fluctuations of density are not important and this 
function can be written as follows: 



(n w (> qthr)) 



1 



3 erfefafer) (l + z) 2 
8 l v (q w (qthr)) ' 

4y/nr)thr eii(rjthr) + 2 exp(-77 t 2 /lr ) 
7rexp(77 2 , lI ,)erfc(77i/ l ,.)[l + erffohr)] 



(68) 



where r] 2 hr — qthr/8r 2 , and the factor (1 + z) 2 describes 
the expansion of the universe. For qthr — > 0, 77^,, < 1, we 



M. Demiariski, and A.G. Doroshkevich: Statistical characteristics of Large Scale Structure 



15 



have (q w (qthr)) — > (?u>(0)) as given by l|62[l and the mean 
linear number density of pancakes increases as 

(n w (> Qthr)) oc (1 + z) 2 t- 2 oc (1 + zf . 

For r/ thr > 1 we have (q w (q t hr)} ~ Qthr- Similar re- 
lations can also be written for the threshold surface 
density of pancakes. Application of these results to 
the Lyman-a clouds observed in a wide range of red- 
shifts (Demiahski, Doroshkevich & Turchaninov 2003; 
Demiahski, Doroshkevich 2003) confirms that this rela- 
tion correctly describes the observed redshift distribution 
of stronger lines in the Ly-a forest. 

7.5. Coagulation approach 

To describe the nonlinear evolution of walls observed in 
deep galaxy surveys we can also use the ID version of the 
coagulation equation (Smoluchowski 1916; Silk & White 
1978) which can be written in the comoving space as fol- 
lows: 



dn(q, t) _ 1 f q 
dr ~ 



dxP(x,q — x,T)n(x,T)n(q — x,t)— (69) 



n(q,r) / dxP(x,q,T)n(x,T) . 
Jo 

Here n(q, r) is the comoving linear number density of walls 
with the dimensionlcss surface density q, P{x, q, r) is the 
rate of aggregation of walls, and it is assumed that walls 
accumulate the main fraction of mass. Using this approach 
it is possible to find the linear number density of walls but 
their evolution depends on the unknown aggregation rate 
P(x, q, t) which is a complicated function of q & r. In 
particular, it depends upon the initial power spectrum. 

The simplest reasonable solution of the coagulation 
equation, similar to l|62|l . can be written as follows: 

n 0?' r ) = , 02^4 ex P ( ~p~2 I > P ( X >Q' T ) = F ot, (70) 



/ P 2 T 4 



and for the mean linear number density of walls we get: 
' 2 ' > ■•' ■ (71) 



{n{qthr, t)) = (1 + z) I dx n(x, r) 
J qthr 



4(1 + zf ( Aq thr 
- l v P r 2 eXP { Pot* 

what is similar to Q68[). For more complicated aggregation 
rates P(x, q, t) some solutions were given by Silk & White 
(1978). 

7.6. Linear number density of low mass pancakes 

The approach discussed in Sec. 7.4 neglects the influence of 
small scale perturbations and approximately characterizes 
only the mean linear number density of richer walls. For 
low mass pancakes of a size comparable with the coher- 
ent scale of initial density field, q w ~ go, the mean linear 



number density of pancakes with a threshold size qthr de- 
pends upon the spectral moment, mo, and go ©■ It can be 
found with the standard technique (see, e.g., BBKS) used 
to describe the condition that a random function exceeds 
a certain value, we get: 

V3(l + z) 2 e(y) 



n{> qthr)) ~ 
erf(»7th r ) 



$ = 



Vthr 



l6Trl v T(z)yfqo~ 
" 2 < r ' 3 dx 2 dx 3 



(72) 



e 2 (y) = i + 



go d 2 £ v _ y 2 
3 d 2 q " 



exp(-T) thr - x: 

A , (p_l)2 



4)- 



1 + 



1 



p(l+p) 2 p 2 (l+p) 



V = qthr/qo, P = y/l + y 2 , 



Vthr = qthr/&T 2 (z) . 



Here the factors (1 + z) 2 and erf(r] t h r ) describe the im- 
pact of expansion of the universe and merging of pan- 
cakes. The factor 0(y) in (f?^ introduces corrections for 
the case q < qo, while 8 — * 1, for q ^> q , y ^> 1. The 
density l|72|l depends upon transverse motions character- 
ized by the parameters rj 2 ,r] 3 . For 772,7/3 > 1, Vthr -C 1, 
(n(> qthr)) oc (1 + z) 3 and it grows not so fast as for 
the richer walls (|68|l because formation of pancakes with 
q < Qo is suppressed. 

The relation l|72|l characterizes pancakes by their 
threshold size at the redshift of formation and neglects 
evolution of pancakes after they are formed. However, the 
surface density of formed pancakes is changing because 
of their transversal compression and/or expansion which 
shifts some of the pancakes under and/or over the obser- 
vational threshold. This problem is quite similar to that 
discussed in Sec. 7.3, where it was noticed that the surface 
density is a more adequate characteristic of DM pancakes 
which takes into account these variations. 

In the relations 1|72JI the threshold size of pancakes, 
qthr, together with their transverse sizes appears only in 
the function <&(r]thr, V21 This means that to go from 
the threshold size to the threshold surface density we must 
link q t hr and <j t hr with the expression l|65|l and find a new 
function $(cr t ^ r , 772, 773). This procedure is quite similar to 
that used in Sec. 7.3, and, for 772 —> 00, i]3 °°, we § e t 
instead of 172|) that 



(n(> (Tthr)) 



V3(i + z) 2 e 

16nl v T(z)^/qo 



exp(-?? 2 ) 



erf(7y s ) 



(73) 



Othr 



8r 2 (l + z) 2 ' 

This result demonstrates again that, as was discussed in 
Sec. 7.3, the transverse compression and/or expansion of 
pancakes compensate each other and it does not change 
their statistical characteristics, if we consider each pan- 
cake as formed at its current redshift. 

Both expressions (|72f) and (|73|l were used in 
Demiahski, Doroshkevich & Turchaninov (2003) to de- 
scribe the observed evolution of the mean linear number 
density of pancakes and to estimate the important param- 
eter qo and the moment mo of the initial power spectrum. 
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8. Summary and discussion 

In this paper we continue the statistical description of 
the process of LSS formation and evolution based on the 
Zel'dovich theory of nonlinear gravitational instability. 
First results obtained in DD99, DDMT and Demiahski 
et al. (2003) show a significant potential of this approach. 
Here we are allowing for deformation of pancakes after 
their collapse along the axis of the most rapid compres- 
sion, interaction of large and small scale perturbations and 
the impact of small scale cutoff in the initial power spec- 
trum. This extension allows one to consider three impor- 
tant problems. 

First of all, we are able to find the mass functions 
and fractions of matter accumulated by the LSS elements, 
namely, pancakes, filaments and halos for a wide range of 
redshifts. As was shown in Sec. 6, the mass functions de- 
scribe reasonably well simulated mass distributions at all 
redshifts that emphasizes the generic character of the pro- 
cesses of formation of all structure elements. These func- 
tions provide quantitative description of the LSS evolution 
that in itself is an important problem. 

Secondly, we discuss the interaction of large and small 
scale perturbations, which manifests itself as a strong con- 
centration of galaxies within filaments and walls observed 
at small redshifts. This interaction is complex and it re- 
quires more detailed investigation. 

Thirdly, we derived the mass function and the mean 
linear number density of pancakes at high redshifts. Both 
functions play an important role in the interpretation of 
the Ly-a forest observed in spectra of the farthest quasars. 
Results obtained in Sec. 7 are successfully applied in 
Demiahski et al. (2003) and Demiahski & Doroshkevich 
(2003) for detailed description of observed absorbers and, 
in particular, lead to estimates of the spectral moment 
mo and the mass of the dominant fraction of dark matter 
particles. 

The rapid growth of the observed concentration of neu- 
tral hydrogen at redshifts z ~ 6 (Djorgovski et al. 2001; 
Becker et al. 2001; Pentericci et al. 2001; Fan et al. 2001) 
is an evidence in favor of the reionization of the Universe 
at this redshift. These observations stimulate discussions 
of the reheating of the universe and, in particular, of the 
warm dark matter (WDM) models (see, e.g., Barkana et 
al. 2001; Loeb & Barkana 2001). This means that it is 
worthwhile to find direct estimates of the small scale ini- 
tial power spectrum and its influence on the LSS forma- 
tion. Results obtained in Sees. 4, 5 & 7 are quite important 
for such investigations. 



8.1. Amplitude of perturbations 

In the Zel'dovich theory the evolution of structure can be 
suitably characterized by an effective dimensionless 'time', 
t(z, Q m ,h), introduced by i|13fl : 



t(z) = T (f2 m , h)B{VL m , z), t 



which describes the evolutionary stage achieved in the 
model. The function B(z), the amplitude A and the spec- 
tral moment m_2 were introduced in 10-01). 

To can be measured by different methods. More popu- 
lar characteristics of the amplitude are <7g - the variance of 
mass in a randomly placed sphere of radius 8/i _1 Mpc, Tq 
- the amplitude of the quadrupole component of the CMB 
anisotropy, and the correlation function of observed or 
simulated matter distribution. Some other methods were 
discussed in DDMT. 

The amplitude of initial perturbations, A, is simply 
linked with ag 



= 9 A 2 



dxx 3 T 2 (x) 



smxg — xs cos £8 



(75) 



where x = k/ko, Xg = 8xfl m h and T(x) was introduced 
in Sec. 2 . Using the latest estimates (Spergel et al. 2003) 
for the ACDM model © we get: 



cr 8 » 0.055A re 0.9 ± 0.1, A S3 16.4 ± 1.82 , 
t = t 8 « (0.22 ± 0.02) . 



(76) 



Using results of Bunn & White (1997), we can link the 
amplitude To with the quadrupole anisotropy of the CMB, 
Tq. For the same basic model we have 



t = t t - 0.27 



0.65 



0.8 



Qrn.h 2 \ Tq 



0.2 J 20fiK 



(77) 



Similar expression can be also written for hyperbolic cos- 
mological models (DD99). 

The variance of displacement (7), a s and t, can be 
directly expressed through the observed two point corre- 
lation function of galaxies, £ ffa i( r )i as follows: 



= lim / dx 



(l - H Xigal[x) . 



and for the correlation function £(r) approximated by the 
power law, 



r < 



n ■ 



we have 



2—7 7 



(2- 7 )(3- 7 ) 



(78) 



Here is the first zero-point of the correlation function. 
It is usually found with a small precision, but for 7 sa 1.5 
- 1.7, 1 — 7/2 S3 0.25 - 0.15, even essential variations of 
do not change significantly the final estimates of a s and 
t. Evidently, the nonlinear clustering of galaxies at small 
r increases the estimate (|78fl of a s . However, analysis of 
simulations (DFTT) shows that this expression underes- 
timates the amplitude of perturbations. 

The correlation functions for the APM survey were 
found in Loveday et al. (1995), for the Las Campanas 
Redshift Survey in Jing et al. (1998) and for the 
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Durham/UKST Redshift Survey in Ratcliffe et al. (1998). 
Using these results we have for (40±10)/i. -1 Mpc that 

t =t c « (0.2 ±0.04) (j^A • (79) 

Applying the relations (Mi2t to the systems of walls se- 
lected in the LCRS, Durham/UKST Redshift Survey and 
the SDSS EDR (Doroshkevich et al. 2002) allows one to 
estimate To as 

r =t w = (0.27 ±0.04)vUnft 2 /0.2, (80) 

what is consistent with the above estimate. 

Differences between estimates of tt, T£ and t w given 
by 1(771 1791 180(1 reflect the precision actually achieved in 
modern observations and, in particular, indicate a possi- 
ble influence of cosmic variance, nonlincarity and random 
factors. 

8.2. Mass function of structure elements 

The Press-Schechter formalism derived for the spherical 
collapse describes quite well the mass functions of vari- 
ous observed and simulated structure elements. In spite of 
this, all attempts to extend it to the collapse of asymmet- 
ric objects failed with the exception of recently proposed 
description of collapse of elliptical clouds (see, e.g., Loeb 
and Barkana 2001; Sheth and Tormen 2002). In this pa- 
per we demonstrate that the Zel'dovich theory successfully 
describes collapse of any object. 

Indeed, formation and relaxation of pancakes along the 
axis of the most rapid compression does not prevent their 
deformation in transverse directions due to relatively small 
gradients of density and pressure in these directions. The 
growth of density within the pancake changes the rate 
of evolution and accelerates its compression and further 
transformation into high density filaments and clouds. 
However, the masses of filaments and clouds formed due 
to such compression remain almost unchanged. The same 
factor also decelerates the expansion of the pancake and 
its dissipation and, so, increases the fraction of surviving 
pancakes. 

This means that, allowing for these deformations, we 
obtain approximate time dependent mass functions of 
structure elements formed due to successive compression 
in one, two and three directions. These results emphasize 
the generic character of the formation of all structure el- 
ements and link the fundamental characteristics of struc- 
ture with the initial power spectrum. They extend the 
Press-Schechter formalism for all LSS elements including 
the filaments and walls which are far from equilibrium. 
The unexpected similarity of the mass functions for these 
elements verifies that the shape of collapsed clouds influ- 
ences the rate of collapse but it does not change signifi- 
cantly their mass functions. 

As is seen from the comparison of results presented 
in Sec. 5 and in Loeb & Barkana (2001) for the Press- 
Schechter formalism, the mass functions derived in both 



approaches are similar and, in particular, the relations 
(I44|l and (|50|l resemble each other. Both mass functions 
are sensitive to the damping of small scale perturbations 
caused by the random motion of DM particles and the 
Jeans damping. They predict the existence of numerous 
low mass objects which can be identified with isolated 
dwarf galaxies and a rich population of Ly-a absorbers. 
Both approaches predict a strong suppression of forma- 
tion of isolated low mass objects with M < M p . 

This comparison indicates that the differences between 
these approaches are quantitative rather than qualitative. 
The similarity of the mass functions demonstrates, in fact, 
the self similar character of the process of structure evolu- 
tion that is the successive condensation of matter within 
clouds, filaments and walls with progressively increased 
sizes and masses. At the same time, the approximate 
character of both the Zel'dovich theory and the Press- 
Schechter formalism implies that the proposed mass func- 
tions are only approximate. Detailed analysis of high res- 
olution simulations with application of the approach pro- 
posed in this paper will allow to improve results obtained 
in Sec. 5. 

Due to their high overdensity above the mean den- 
sity, both filaments and walls are easily detected in the 
deep galaxy surveys using the Minkowski Functional ap- 
proach (Schmalzing et al. 1999; Kerscher 2000) and the 
well known friend-of-friend method generalized in the 
Minimal Spanning Tree technique (Barrow, Bhavsar & 
Sonoda 1985; van de Weygaert 1991; Doroshkevich et al. 
2001, 2002). The mass functions of the observed struc- 
ture elements can be compared with the expected ones. 
For simulations the same comparison can be performed 
at all redshifts what allows to trace the expected redshift 
dependence of these functions. 

In Sec. 5.7 for the first time we compare the expected 
mass functions with simulated ones. In Doroshkevich et al. 
(2002) the expected mass functions are compared with the 
observed ones for the LSS elements selected with various 
threshold overdensity from the SDSS EDR. These results 
show that both relations (|4*5|) and (|5U|l describe quite well 
even the mass distribution of filaments and walls with a 
moderate richness. Stronger disagreement appears for the 
richest walls and filaments formed due to the process of 
percolation which is not described by the Zel'dovich the- 
ory. 

However, the potential of both approaches is limited 
as they cannot describe the final disruption of collapsed 
clouds what leads to formation of numerous low mass 
satellites of the central object and incorporation of fila- 
ments and walls into a joint network (percolation process). 
Neither the Press-Schechter formalism nor the Zel'dovich 
theory can describe the impact of environment discussed 
in Sec. 4 and, in particular, the faster evolution of clouds 
accumulated by richer walls as compared with isolated 
clouds. As was discussed in Sec. 5, both approaches do 
not describe the coagulation of high density clouds what 
distorts their mass function and makes it similar to the 
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power low. These problems remain open for further inves- 
tigations. 

8.3. Interaction of large and small scale perturbations. 
Large scale bias 

The analysis of redshift distribution of absorbers shows 
almost homogeneous spatial distribution of baryonic and 
DM components of the matter on scales > lft. _1 Mpc. At 
the same time, the observed spatial distribution of lu- 
minous matter is strongly nonhomogeneous. Thus, about 
half of the observed galaxies are concentrated within large 
walls while within the Bootes void the number of galaxies 
is very small. The observed walls and voids are associated 
with compressed and expanded regions of the Universe 
and these observations point out to the possible correla- 
tions between the galaxy formation and the large scale 
deformation field (see, e.g., Rees 1985; Dekel & Silk 1986; 
Dekel & Rees 1987). Some of such correlations were al- 
ready noticed in simulations (see, e.g., Sahni et al. 1994). 

As was shown in Sec. 4, the Zel'dovich theory allows 
to quantify this interaction and demonstrates that the for- 
mation rate of high density objects is modulated by large 
scale perturbations. This modulation cannot significantly 
change the redshift evolution of the fraction of mass accu- 
mulated by structure elements. However, this interaction 
accelerates the formation of high density halos and galax- 
ies within regions associated with the future LSS elements 
such as clusters of galaxies, filaments and walls. These re- 
sults suggest that the poorer sample of isolated galaxies 
and invisible DM halos situated within low density regions 
were formed later then similar galaxies and DM halos sit- 
uated within filaments and walls. The results presented 
in Sec. 4.1 illustrate this influence for simulated matter 
distribution. 

Another important factor which can not be neglected 
is the acceleration of merging and coagulation processes 
within filaments and walls. Both processes change the 
mass function of observed galaxies and increase the frac- 
tion of massive objects. The high concentration of ob- 
served galaxies within filaments and walls and small frac- 
tion of isolated galaxies (< 10%) is an evidence in favor 
of the large scale bias. These factors also determine the 
observed correlation function of galaxies. 

8.4. Properties of low mass pancakes and the mass of 
DM particles 

Numerous simulations performed recently (see, e.g., Zhang 
et al. 1998;, Weinberg et al. 1999) demonstrate that the 
absorption lines observed in spectra of the farthest quasars 
- the Ly-tv forest - are related to the numerous low 
mass clouds formed at high redshifts. As was discussed in 
Demianski et al. (2003), some of the statistical character- 
istics of observed absorbers can be successfully explained 
on the basis of Zel'dovich theory. At the same time, the 
analysis shows that the approximate description based on 



results obtained in DD99 have to be essentially improved. 
Indeed, the absorbers are observed in a wide range of red- 
shifts and the properties of long-lived absorbers are chang- 
ing with time at least due to their transverse expansion 
and compression. 

Discussion of this problem in Sec. 7 shows that for sta- 
tistically homogeneous sample of absorbers we can neglect 
this influence, at least for their two important characteris- 
tics, namely, for the PDF of their surface density l(f)2j) and 
their mean linear number density (|72|l . Results obtained 
in Sec. 7 link these characteristics with the properties of 
the initial power spectrum at small scale. This means that 
the analysis of these characteristics allows one to measure 
the variance of initial density perturbations and to restrict 
the mass of dominant fraction of dark matter particles. 

Results obtained in Sec. 7 we applied to the Lyman-a 
forest in Demianski, Doroshkevich & Turchaninov (2003). 
This analysis confirms that the relations (|62|l and l|72|l 
describe quite well the PDF and the redshift distribution 
of ~ 5000 observed absorbers. These results also verify 
the Gaussianity of initial perturbations. In turn, analysis 
of redshift distribution of absorbers allows one to estimate 
the actual characteristics of the initial power spectrum as 
follows: 

q w 0.01 ± 0.003, m « 0.2 - 0.5, M DM > 1 - 5keV(81) 

Analysis of simulations (Narayanan et al. 2000; Barkana, 
Haiman & Ostriker 2001) restricts the mass of DM parti- 
cles to Mdm > 1-1.5 keV. Comparison of other character- 
istics of pancakes derived in Sees. 3 & 7 with observations 
can be found in Demianski, Doroshkevich & Turchaninov 
(2003). 

8.5. Reheating of the universe 

The relations (|28l 142(1 show that in the Zel'dovich theory 
of gravitational instability the rate of matter collapse at 
high redshifts depends upon both the amplitude of pertur- 
bations, t(z), as given by (|13|) . and the coherent length of 
density field, l p & qo, as given by which, in turn, de- 
pends on the shape of the initial power spectrum at large 
k and the mass of the dominant type of DM particles. 

As was noticed in Sec. 3, the most effective matter 
condensation within high density clouds occurs at r(z r ) = 
r r w (see equation (J22J) when these clouds already 

accumulate ~ 2% of matter and the main fraction of mass 
is concentrated within clouds with M ~ (0.2 — 0.5) (M). 

For go ~ 10 -2 we can estimate the redshift of the pe- 
riod of most efficient condensation and the mean mass of 
DM clouds as follows: 

T r ~ 0.04, 1 + z r w 1.3to/tv - 7- 8, 

(M DM ) ~ r?M v « 1O 12 M Q^Y f^V , (82) 

what is similar to the mass of a typical galaxy. For qo ~ 
10 -3 we have, respectively, 

T r - 0.013, 1 + z r « 1.3t /tv - 20 , 
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(M DM ) ~ 4 • 10~ 2 t?M v « 1O 9 M f ^) - (83) 

what is similar to the mass of a typical dwarf galaxy. These 
results illustrate strong links between the period of reheat- 
ing and the shape of the power spectrum at small scale 
characterized by qo- Even at higher redshifts noticeable 
fraction of matter can be compressed within high density 
clouds with masses of galactic scales. 

According to available estimates (Haiman & Loeb 
1999; Loeb and Barkana 2001) the reionization of the uni- 
verse occurs after ~ 1 - 3% of matter is concentrated 
within high density halos. From (|82l I83|) it follows that 
this fraction of collapsed matter can be reached already 
at r ~ r r at redshifts z = z r ~ 8 - 25, for qo ~ 10 -2 — 10~ 3 . 
The effective reionization and reheating of the universe at 
such redshifts is consistent with the observed concentra- 
tion of neutral hydrogen at redshifts z ~ 6, which charac- 
terizes mainly the rate of generation and achieved inten- 
sity of UV background. This means that these ranges of 
qo and z r provide a reasonable estimate of the period of 
reheating. 

The observations of environment of high redshift 
quasars (Fan et al. 2001) provide an evidence in favor of re- 
heating of the Universe at z ~ 6 - 10, what is more consis- 
tent with qo r~j 10~ 2 . However, our analysis shows that, for 
the standard WDM model with the Harrison - Zel'dovich 
large scale power spectrum and exponential cutoff caused 
by a mass Moai ~ 1 kcV of the dominant fraction of DM 
particles, some problems can appear with formation of low 
mass isolated galaxies with M < 10 8 — 10 9 M Q . Perhaps, 
some excess of power at small scales can help to solve this 
problem. First observational indications of such excess are 
presented in Demiahski & Doroshkevich (2003) 
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Appendix A: Correlation functions of initial 
perturbations. 



Appendix B: Characteristics of the deformation 
field. 



In this appendix we present a few correlation and structure The general deformation of a spherical cloud with a diam- 

functions which describe the relative spatial distribution of eter q = q/l v can be suitably characterized by the dimen- 

important parameters of the initial perturbations. These sionless random scalar function 
functions have been introduced in DD99, where more de- 
tails can be found. 



e(q) = 0(g,M) = 



\/3[S(q/2) - S(-q/2)] ■ q 



(T s q 



(B.l) 



The structure function of gravitational potential per- 
i ,. i , • , ,. r , , ., ,. , instead of the deformation tensor ct,t. Expansion of this 

turbations characterizes correlation or the gravitational r 



potential in two points qi and q 2 . As the power spectrum 



function into spherical harmonics characterizes deforma- 



c , . fiiiuij- i r u tions of the cloud with a required accuracy. This problem 

is a tunction of only the absolute value of wave number ... , . . 

1,1,,.; , f ,. , j - i~ ~i i is similar to the usual description of CMB anisotropy in 

|fc|, this structure function depends on gi 2 = |qi — q 2 | and ... 



for the perturbations of gravitational potential we have 
< A0i 2 A0 34 > 



l 2 a 2 



G a {qu) - Go (<?i3)+ 



(A.l) 



G (q 23 ) - GofeO, 9ij = |q< - Oil i 

A0i2 = 0(qi) - 0(q 2 ), A034 = 0(q 3 ) - 0(q 4 ) 



Go(q) 



p(k) 



sin kq 
kq 



dk. 



Differentiations of this structure function give other struc- 
ture functions: 

3 < A0(q 12 )g,(q3) > = 
l v a 2 

(qi - q3)iGi(<zi 3 ) - (q2 - qs^Gife) > 

^ = SijGifau) + (qi2Mqi2)jG 2 (<7i 2 ) , 



Gi(a;) = G /z, G 2 (x)=G' 1 /x, G 3 (x)=G' 2 /x, (A.2) 

G 12 (z) = Gr(x) + a; 2 G 2 (x) = G ' = (xGj)' , 

For the Harrison - Zel'dovich primordial power spec- 
trum and for the CDM transfer function (BBKS), these 
functions can be approximated by 



Gi(x) m (1 + q )[l + yj 'q 2 + x 2 + a x 2 ] \ 
G 12 w G?(l + go)" 1 !! - *0X 2 + qll^ql+x 2 ] , 

Gi(0) = G 12 (0) = 1, G 2 (0) = -[g (l + go)]" 1 . 
x = q/l v , qo = l p /l v , a a = 5 (l v /L ) 2 » 0.3 , 
where the typical scale -Lo is defined as 



(A.3) 



Ln = 3tt 2 



P(A) 1 - 



sin fcLo 
kLn 



k~ 2 dk, 



p(k)dk , 



and the scales /„&go were introduced by (TJJ. 



terms of spherical harmonics, and we have: 

0( 9 ,M)=X; Olm(«)nm(e,0), 



(B.2) 
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l,m \ 



= {-i) m i 



l 21 + 1(1- m)\ 



P l m (cos9)e 



irrwj) 



Y, 



(-1)1-™^ 



4-7T (/ + m)!' 

,<f>), Z > to > 0. 



(e( qi )9(q 2 )) = Ml2 ) = a 2 (g) £ & 2 (g)P*(Mi 

I 

Gifo) + Gi(a; 2 ) - Gi 2 (xi) - Gi 2 (x 2 )+ 
Hi2{Gi{x 2 ) - Gx{xx) + G 12 {xi) - Gi 2 (x 2 )] 

qiq2 



A*12 



9192 



•i'i = g 



1 



A* 12 



21- 



x 2 



1 



A*12 



-/ 



4tt 



47T 



*(g,Mi2 = 1) = a 2 (g), *(g,Mi2 = 0) = <r 2 (g)r s (g), (B.4) 

„2 



<7 2 (g)=2[l-G 12 (g)], r s (g) = -\G 2 



\V2 



(B.5) 



Here cr 2 (q) and r s (g) are the dispersion and the coeffi- 
cient of correlation of dimensionless orthogonal displace- 
ments, and the functions Gi, G 2 and Gi 2 were introduced 
in Appendix A. The correlation coefficient r s is a weak 
function of q and for q -C 1, go <C 1, we have 

1/3 < r s < 2~ 3/2 , (B.6) 

what is consistent with the usual correlations of the com- 
ponents of deformation tensor. 

The function ^(q,fii 2 ) is symmetrical with respect to 



the replacement \i\ 2 



-/ii2 and, therefore, bf — for 



odd I. Using the relations (|A.3(1 we have for the relative 
amplitude of even spherical harmonics, for the most inter- 
esting cases qo <C q -Cl and q <C <?o -Cl, respectively: 



0.533, 6, w 0.381, 



& 2 « 0.037, 6g 



0.014, ... (B.7) 
(B.8) 

These results justify the assumption used in DD99 to ne- 
glect higher order harmonics of perturbations with I > 4. 



^5 



0.55, 



b% 



0.44 > bj,b 2 . 



